xi_func=@(rho)(1-rq_ratio)*max([1+rho,1/(1+pi_m)])+rq_ratio/(1+pi_m)-cost;
foc_bank=@(z,rho)xi_func(rho)*(z*ddinv(z)/Nb+dinv(z))-1;
rho_bar = fzero(@(x)xi_func(x)-1/beta,[1e-10,2]);
rho_grid=linspace(1/(1+pi_m)-1,rho_bar,100);

func_var.beta=beta;
func_var.alpha1=alpha1;
func_var.alpha2=alpha2;
func_var.alpha3=alpha3;
func_var.lambda_func=lambda_func;
func_var.dlambda_func=dlambda_func;
func_var.rq_ratio=rq_ratio;
func_var.cost=cost;
func_var.p_bar=p_grid(end);
func_var.rho_bar = rho_bar;
%func_var.Nb=Nb;